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We review the theoretical status of squark and gluino hadroproduction and provide nu- 
merical predictions for all squark and gluino pair-production processes at the Tevatron 
and at the LHC, with a particular emphasis on proton-proton collisions at 7 TeV. Our 
predictions include next-to-leading order supersymmetric QCD corrections and the re- 
summation of soft gluon emission at next-to-leading-logarithmic accuracy. We discuss 
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the impact of the higher-order corrections on total cross sections, and provide an esti- 
mate of the theoretical uncertainty due to scale variation and the parton distribution 
functions. 
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1. Introduction 

The search for supersymmetry (SUSYjl^ is a central part of the experimental 
program at the hadron coUiders Tevatron and LHC. Models of weak-scale SUSY 
provide a promising solution to the hierarchy problem of the Standard Model (SM) 
and comprise new supersymmetric particles (sparticles) with masses of order 1 TeV. 
The coloured TeV-scale sparticles, squarks {q) and gluinos [g), would be produced 
copiously in hadronic collisions and thus offer the strongest sensitivity for super- 
symmetry searches at the Tevatron and the LHC. 

We cons ider the minimal supersymmetric extension of the Standard Model 
(MSSMj2MI where, as a consequence of R-parity conservation, squarks and gluinos 
are pair-produced in collisions of two hadrons hi and /12: 

hih2 qq ,qq ,qg ,gg, tiii, 1212 + X . (1) 

The production of top squarks (stops), ^1.2, has to be treated separately, because 
the strong Yukawa coupling between top quarks, stops and Higgs fields gives rise 
to potentially large mixing effects and mass splitting.^ In Eq. ([1]) and throughout 
the rest of this paper, ti and t2 denote the lighter and heavier stop mass eigenstate, 
respectively. For the other squarks we suppress the chiralities, i.e. q = iqL,qR)-, and 
do not explicitly state the charge-conjugated processes. 

Searches for squarks and gluinos at the proton-antiproton collider Tevatron with 
a centre-of-mass energy of ^/S = 1.96 TeV have placed lower limits on squark and 
gluino masse s in the range of 300 to 400 GeV, depending in detail on the specific 
SUSY model.l^^^ The proton-proton collider LHC, which has been operating at 
^/S = 7 TeV in 2010, has already significantly extended the squark and gluino mass 
limits to values of aro und 85 GeV.I^HISl Dedicated searches for the lighter stop 
mass eigenstate at and the TevatrorP^m^ have placed lower limits in 

the range 70 to 200 GeV. Already in 2011, with a projected integrated luminosity of 
1 to 2fb~^. the LHC should be sensitive to squarks and gluinos with masses in the 
TeV region!^, while SUSY particles with masses u p to 3 TeV can be probed once 
the LHC reaches its design energy oi^=U TeV.lSQiail 

Accurate theoretical predictions for inclusive squark and gluino cross sections 
are needed both to set exclusion limits and, iri case SUSY is discovered, to determine 
SUSY particle masses and properties .^2211251 ^pj-^g inclusion of higher-order SUSY- 
QCD corrections significantly reduces the renormalization- and factorization-scale 
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dependence of the predictions. In general, the corrections a lso incr ease the size of 
the cross section with respect to the leading-order predictiorPSHU! jf the renormal- 
ization and factorization scales are chosen close to the average mass of the produced 
SUSY particles. Consequently, the SUSY-QCD corrections have a substantial im- 
pact on the determination of mass exclusion limits and would lead to a significant 
reduction of errors on SUSY mass or parameter values in the case of discovery. 
The processes listed in Eq. ([T]) have been known for quite some time at next-to- 
leading order (NLOlin SUSY-QCD .EHHSS] Electroweak corrections to the ©(a^) 
tree-level proces sea^^H^ and the electroweak Born production channels of 0{aas) 
and ©(a^pnUH g^]-g general significant for the pair production of SU(2)-doublet 
squarks and at large invariant masses, but they are moderate for inclusive cross 
sections. 

A significant part of the NLO QCD corrections can be attributed to the threshold 
region, where the partonic centre-of-mass energy is close to the kinematic production 
threshold. In this region the NLO corrections are dominated by soft gluon emission 
off the coloured particles in the initial and final state and by the Coulomb corrections 
due to the exchange of gluons between the massive sparticles in the final state. The 
soft-gluon corrections can be taken into account to all orders in perturbation theory 
by means of threshold resummation techniques) ■ ' 

Recently, such a threshold resummation has been performed for all MSSM 
squark and gluino production processes, Eq. ([T]), at next-to- leading-logarithmic 
(NLL) accuracy.'^^H^ A formalism has been developed in the framework of ef- 
fective field theories which allows for the resummation of soft and Coulomb gluons 
in the production of c oloured sparticles, but has so far only been applied to squark- 
antisquark production.l^SEHI in addition, the dominant next-to-next-to-leading or- 
der (NNLO) corrections, including those coming from the resummed cross section at 
next-to-next-to-leading-logaxithmic (NNLL) level, have been calculated for squark- 
antisquark pair-production.ISniSIl xhe production of gluino bound states as well as 
bound-state effects in gluino-pair production has also been studicd.l^^ES] 

In this work we will present the state-of-the-art SUSY-QCD predictions for 
the MSSM squark and gluino hadroproduction processes, Eq. ([T]), at the Tevatron 
and the LHC, including NLO corrections and NLL threshold resummation. The 
processes will all be treated on the same footing, i.e. classes of beyond-NLO effects 
that have been calculated only for one particular process or reaction channel will 
not be taken into account. We will discuss the impact of the SUSY-QCD corrections 
on the total cross sections and provide an estimate of the theoretical uncertainty 
due to scale variation, parton distribution functions, and the strong coupling as- 

The structure of the paper is as follows. In section [2] we briefiy review the ap- 
plication of the resummation technique to total cross sections for coloured sparticle 
pair-production. The numerical results are presented in section [3l We show predic- 
tions for the Tevatron, and for the LHC with centre-of-mass energies of = 7 TeV 
and y/S =U TeV. We will conclude in section m 
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2. NLL resummation 

In this section we provide a brief background to the calculation of the threshold- 
resummed cross sectionJ33H3Sl -^g ^^se for our results in the next section. The re- 
summation for 2 — > 2 processes with all four external legs carryi ng colo ur has 
been studied extensively in the literature, specifically for heavy-quarli^2l55] ^^^^ jg-j- 
production.'S^'^In our calculations we make use of the framework developed there. 

The hadronic threshold for the inclusive production of two final-state particles 
k, I with masses mk and mi corresponds to a hadronic centre-of-mass energy squared 
that is equal to S' = {rrik -|-m;)^. Thus we define the threshold variable p, measuring 
the distance from threshold in terms of energy fraction, as 

p = ("^^■ + "^')\ (2) 

Our results are based on the following expression for the NLL-resummed cross 

on _ oo 

section, matched to the exact NLO calculation^ 

(NLO+NLL)/ r 21 2\ (NLO) / r 21 2\ 



- E 



dN /,/^, {N + 1, p^) {N + 1, p^) 

CT 



-':;lAN,{m^},p') - al;l,(7V,{mn,A^^) |<.,o,J > (3) 

where the last term in the square brackets denotes the NLL resummed expression 
expanded to NLO. The initial state hadrons are denoted generically as hi and /12, 
and p is the common factorization and renormalization scale. The resummation is 
performed after taking a Mellin transform (indicated by a tilde) of the cross section, 

dhih2'^ki{N,{m^}) = I dp p^^^ ahih2'^ki{p,{'^^}) ■ (4) 
Jo 

To evaluate the contour CT of the inverse Mellin transform in Eq. ([3]) we adopt 
the so-called "minimal prescription" The NLL resummed cross section in Eq. ([3]) 
reads 

X A,(7V + l,g2,M')A,(7V + l,Q2,A*')Ag^,;^,(7V + l,Q2,^2^, (5) 

where o'^^l^/.; / are the colour-decomposed leading-order cross sections in Mellin- 
moment space, with I labelling the possible colour structures! 45|46 | jjgj-g have 
introduced the hard scale = (m^ -I- to;)^. The perturbative functions Cij^kij 
contain information about hard contributions beyond leading order. This informa- 
tion is only relevant beyond NLL accuracy and therefore we keep Cij^kij = 1 in 
our calculations. The functions and Aj sum the effects of the (soft-)collinear 
radiation from the incoming partons. They are process-independent and do not de- 
pend on the colour structures. They contain the leading logarithmic dependence, as 
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well as part of the subleading logarithmic behaviour. The expressions for A.; and Aj 
can be found in the literature!^ The resummation of the soft-gluon contributions, 
which does depend on the colour structures in which the final state SUSY particle 
pairs can be produced, contributes at the NLL level and is summarized by the factor 



Af{N,Q^fi') = exp 



'5/"' dq aM ^ 

Uj 

q TT 



(6) 



The one-loop coefficients Dj follow from the threshold limit of the one-loop soft 
anomalous-dimension matrix P^^^^' 

Two remarks can be made regarding squark-gluino production, and stop-antistop 
production in particular. The former reaction is the only one in our set where heavy, 
coloured SUSY particles of different mass are produced. Our resummed expressions 
are sensitive to these differences throu gh t he Born cross sections, and through the 
resummed exponents at the NLL level.^^ In stop-antistop production through qq 
annihilation the gluino exchange diagram, which would require top parton distri- 
bution functions, is missing in the five-flavour scheme we adopt. As a consequence, 

the Born cross section for stop production is proportional to ^^1 — 4to~/s ^ , as 



opposed to y 1 — 4m?/s for production of the other squark fiavours. Here, rriq is 
a generic squar k m ass and s denotes the partonic centre-of-mass energy squared. 
We have argued^ that this has no effect on the resummed expression other than 
through the explicit expression for the Born function. 



3. Numerical Results 

We present numerical predictions for squark and gluino production at the Tevatron 
{VS = 1.96 TeV) and the LHC {VS = 7 and 14 TeV). We compare LO, NLO 
and NLO+NLL matched results, and discuss the theoretical uncertainty due to the 
choice of renormalization and factorization scales and due to the uncertainty in 
the parton distributio n fu nctions (pdfs) and the QCD coupling as- The LO and 
NLO cross sectionJ^SHSSl ^re available in the form of the public computer code 
Prospino.l^ The MS-scheme with five active flavours is used to define as and the 
parton distribution functions at NLO. The masses of the squarks and gluinos are 
renormalized in the on-shell scheme, and the SUSY particles are decoupled from 
the running of as and the pdfs. 

As mentioned in the introduction, the production of stops has to be treated 
separately because of potentially large mixing effects and mass splitting. The pro- 
duction of the other squark flavours, which we assume to be mass degenerate, is 
treated together, i.e. we sum over five fiavours of squarks, q £ {u, d, c, s, 6}. In 
that case our numerical predictions include both chiralities (qL and qn) and the 
charge- conjugated processes. 

The stop cross section is shown separately. Note that since mixing in the stop 
sector enters explicitly only through higher-order diagrams, the stop-mixing angle 
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9f need not be renormalized and one can use the lowest-order expression derived 
from the stop mass matrix. Beyond LO the stop cross section does not only depend 
on the stop mass, but also on the gluino mass rrig, the average mass of the first 
and second generation squarks rriq and the mixing angle 6f. We have fixed rrig, 
VTLq and Of according to the SPSla' benchmark scenario!^ Note, however, that the 
dependence of the stop cross section on the SUSY parameters that enter only at 
NLO is numerically very small with variations of at most 2%!33 The numerical 
results presented for stop production also apply to sbottom production when the 
same input parameters are adopted, since the impac t o f bottom-quark induced 
contributions to sbottom hadroproduction is negligible!^ 

For convenience we define the average mass of the final-state sparticle pair m = 
(m/c -I- mi)/2, which reduces to the squark and gluino mass for qq, qq, and gg 
final states, respectively. The renormalization and factorization scales are taken 
to be equal, /ir — fip ~ fi. As our default, hadronic NLO and NLO-I-NLL cross 
sections are obtained with the 2008 NLO MSTW pdfs" and the corresponding 
as(Mz) = 0.120. In particular, aU plots show resuhs for the MSTW pdf set.^^ We 
also present resuhs based on the CTldS^I and CTEQ6LlE5l pdf sets in tables, for 
comparison. 

Let us now discuss the numerical results for squark and gluino hadroproduction 
at the Tevatron, and at the LHC operating with 7 and 14TeV hadronic centre-of- 
mass energy. We shall study the scale dependence of the LO, NLO and NLO+NLL 
cross sections, the impact of the NLL threshold resummation, and present our best 
predictions at NLO-fNLL for the inclusive cross sections, including the theoretical 
uncertainties from scale variation as well as the pdf and as errors. We put special 
emphasis on the predictions for the LHC at 7 TeV energy, which are of immediate 
importance for the upcoming SUSY searches, and collect detailed results for rep- 
resentative LO, NLO and NLO-I-NLL cross sections and the corresponding theory 
uncertainties in tables. 



3.1. Tevatron 

The state-of-the-art NLO-I-NLL SUSY-QCD cross-section predictions for the indi- 
vidual processes listed in Eq. ([1]), occuring inpp collisions at the Tevatron, are shown 
in Fig.[T] as a function of the average mass m of the final state sparticles. For illus- 
tration we show these results for the case rriq = rrig. In Fig. [2] the total NLO-I-NLL 
cross section for the sum of all four processes, i.e. qq ,qq , qg and gg production, is 
presented. For the central values, the renormalization and factorization scales are 
taken as /i = /ig = m, which is the scale choice a dopt ed as the preferred one in the 
context of the NLO SUSY-QCD calculations .I^SHMl The error band represents an 
estimate of the theoretical uncertainty of the cross section prediction. It consists 
of the 68% C.L. pdf and ctg error, added in quadrature, and the error from scale 
variation in the range m/2 < /i < 2m added linearly to the combined pdf and as 
uncertainty. By this linear combination of scale uncertainty and combined pdf and 



January 20, 2013 14:54 



susy-review 



Squark and gluino hadroproduction 7 

as errors we provide a conservative estimate of the theory error. 

At the scale fj, = m the cross-section predictions are enhanced by soft-gluon 
resummation. The relative X-factor -STnll = cnlo+nll/cnlo for this scale choice 
is displayed in Fig.|3] for squark and gluino masses in the range between 200 and 
600 GeV, and stop masses between 100 and 300 GeV. The soft-gluon corrections 
are moderate for qq and titi production, but increase the predictions for qq, gg 
and qg final states by around 15, 20 and 40%, respectively, assuming squark and 
gluino masses near 500 GeV. Because of the increasing importance of the threshold 
region, the corrections in general become larger for increasing sparticle masses. 
The large effect of soft-gluon resummation for qg and gg production can be mostly 
attributed to the importance of gluon initial states for these processes. Furthermore, 
the presence of gluinos in the final state results in an enhancement of the NLL 
contributions, since in this case the Casimir invariants that enter the Dj coefhcients 
in Eq. (j6|) are larger than for processes involving only squarks. The substantial value 
of -R'nll for qq production at the Tevatron is a consequence of the behaviour of the 
corresponding NLO correctional, which strongly decrease with increasing squark 
mass. 

Let us next discuss the scale dependence of the SUSY- Q CD cross-section pre- 
diction in some more detail. Fig.|4] shows the scale dependence in LO, NLO and 
NLO-I-NLL for the different production processes, listed in Eq. ([Ij, at the Teva- 
tron. The mass of the (w, d, c, s, 6)-squarks and the gluino mass have been set to 
rriq = rrig = 500 GeV, while the stop mass is fixed to rrif^ = 200 GeV. Note that the 
LO predictions are obtained with the LO MSTW pdf selP^ and the corresponding 
LO value for as- The renormalization and factorization scales are set equal to each 
other and varied around the average mass of the final state sparticles, fiQ = m. We 
observe the usual strong reduction of the scale dependence when going from LO 
to NLO. A further significant improvement is obtained when the resummation of 
threshold logarithms is included, in particular for gg and qg production. For the 
qq , qg and gg final states, contributing the most to the inclusive squark and gluino 
cross section at the Tevatron, the LO, NLO and NLO-I-NLL cross section predic- 
tions converge particularly well near = m/2, which appears to be the preferred 
scale choice for these processes. 

3.2. LHC @ 7TeV 

SUSY searches at the LHC, which is currently operating at 7TeV, will soon be 
sensitive to sparticles with masses in the TeV-range. We thus discuss the cross 
section predictions for the LHC @ 7 TeV in more detail, including in particular the 
dependence of the NLL corrections on the ratio of the gluino and squark masses, 
and the size of the theory uncertainty due to scale, pdf and as errors. 

As before, we first present the NLO-)-NLL cross section predictions for the five 
individual processes and for the sum of the qq , qq , qg , gg final states, cf. Figs. [5] and 
|6l We include the estimate of the theory uncertainty obtained from adding linearly 
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the scale dependence in the range to/2 < ^j, < 2m to the combined 68% C.L. pdf 
and as error, added in quadrature. 

In Fig. [7] we show the enhancement of the cross section due to the NLL resum- 
mation at the scale ^ = m. For gluino-pair and squark-gluino production processes 
we find a significant increase of 10-20% at masses around 1 TeV. An enhancement 
of up to 10% can also be observed for the production of heavy stop particles. Note 
that the singularities at the stop-decay threshold = mt + rrig = 782.5 GeV origi- 
nate from the stop wave-function renormalization. They are an unphysical artefact 
of an on-shell scheme approaclPSl and could be removed by taking into account the 
finite widths of the unstable stops. The effect of the NLL resummation on the cross 
section for inclusive squark and gluino production is shown in Fig.|8l 

Next, we present the scale dependence in LO, NLO and NLO-I-NLL for the 
various production processes, see Fig.[9l Here, the mass of the (u, d, c, s, 5)-squarks 
and the gluino mass have been set to rriq — nig — 700 GeV, while the stop mass is 
fixed to TO(^ = 500 GeV. As anticipated, we observe a significant reduction of the 
scale dependence when going from LO to NLO and from NLO to NLO-I-NLL. Note 
that also for squark and gluino production at the LHC, the convergence of the LO, 
NLO and NLO-fNLL cross section predictions near the scale /i = to/2 is striking. 

As the cross section predictions for the LHC operating at 7 TeV are of partic- 
ular phenomenological importance, let us discuss the results in some more detail. 
First, we collect representative values for LO, NLO and NLO+NLL cross sections 
in TablesHEl We include the scale dependence in the range m/2 < fi < 2m for 
the LO, NLO and NLO-fNLL calculations, and the 68% C.L. pdf and as error 
at NLO. The NLO and NLL i^-factors are displayed for convenience. The NLO 
and NLO-I-NLL theory predictions with the default MSTW pdf setP^ are compared 
to those obtained with the CTIO pdf self^. The LO cross sections are calculated 
with the LO MSTWSland the CTEQeLl^^ pdfs. While the LO cross sections dif- 
fer significantly between the two pdf sets, the central NLO+NLL predictions are 
consistent within the theoretical uncertainty. Note, however, that the 68% C.L. pdf 
error estimate for the CTIO set, obtained through rescaling of the 90% CL error 
estimat is significantly larger than that of MSTW. 

The size of the cross section and the impact of the higher-order corrections 
depends on the ratio r = mg/mq. In Fig. [10] we thus present the NLL if-factor for 
different values of r for individual processes with the qq , qq , qg , gg final states . In 
general, in the range of values of r studied here, the dependence of the soft-gluon 
enhancement on r is moderate. Of particular interest is the limit r — >■ 0, i.e. the 
limit of very heavy sq uar ks. Such scenarios are predicted within so-called models 
of split supersymmetr}E3 and lea d to int eresting new phenomenological signatures 
at the LHC, like long-lived gluinos.'^^'^In Fig.[TT]we present the NLO-I-NLL cross 



^The pdf errors presented in the tables apply to individual processes only. In general, the pdf error 
for inclusive squark and gluino production will differ from the one obtained through combination 
of errors for individual processes due to possible correlations between various production channels. 
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section prediction for gluino-pair production in the limit r = and show the impact 
of soft-gluon resummation. The gg cross section for r = is higher than for r — 1, 
with the difference growing from a factor of 1.2 at m = 200 GeV up to a factor of 
4.7 for m = 1.2 TeV. The impact of soft-gluon resummation is significant, and the 
NLL corrections enhance the cross section by about 12% for gluinos with masses 
around 1 TeV. 

The enhancement of the cross section prediction due to higher-order corrections 
depends, of course, on the choice of scale. Choosing /x = m/2, we find in general 
NLL iiT-factors close to one for sparticle cross sections at the LHC with 7 TeV. 
However, a crucial improvement when going from LO to NLO and from NLO to 
NLO-I-NLL is the reduction of the scale dependence. Fig.[T2] shows a comparison 
of the NLO and the NLO-I-NLL scale dependence for the sum of the qq,qq ,qg , gg 
final states as a function of the average sparticle mass, assuming r = 1. Threshold 
resummation leads to a significant reduction of the scale dependence over the full 
range of sparticle masses, with an overall scale uncertainty at NLO-I-NLL of less 
than 10%. In a phenomenological analysis, the different production channels may 
contribute with different weight, depending in detail on the signature and the choice 
of selection cuts. Thus, we plot the scale uncertainty at NLO-I-NLL for the different 
channels in Fig. [131 We also show the full theory uncertainty, consisting of the 
68% C.L. pdf and as error added in quadrature, combined linearly with the scale 
variation error for the NLO cross sections and for the NLO-I-NLL cross sections. We 
find that even though the pdf uncertainty is significant, the inclusion of threshold 
resummation leads to a sizeable reduction of the overall theory uncertainty. This 
is particularly true for the case of gluino-pair and squark-gluino production. For 
gluino-pair production, the total theory uncertainty can be reduced by as much 
as a factor of two when going from NLO to NLO-I-NLL. Looking at all different 
production processes, the overall theory uncertainty at NLO-I-NLL is approximately 
20% or smaller. 

3.3. LHC @ 14 TeV 

SUSY particles with masses in the multi-TeV region can be probed at the LHC 
running at or near its design energy of 14 TeV.^^OET] complete this review and 
to show how the impact of the higher-order corrections depends on the collider 
energy, we also present predictions for squark and gluino production at the LHC 
with 14 TeV. Our cross section predictions at NLO-I-NLL including the full theory 
uncertainty, consisting of the 68% C.L. pdf and as error, added in quadrature, 
combined linearly with the scale variation error, are shown in Figs. [14] and [151 As 
for the LHC at 7 TeV, the impact of NLL threshold resummation is particularly 
significant for gluino-pair and squark-gluino production, with NLL corrections of 
about 30% for gg final states and squark and gluino masses of 2.5 TeV, see Fig.lTSl 
Finally, in Fig.[T7lwe show the LO, NLO and NLO-f NLL scale dependence for squark 
and gluino masses of 1 TeV and ^ = 500 GeV. The results are qualitatively very 
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similar to those obtained for the LHC with 7TeV, i.e. we find a strong reduction 
of the scale dependence when including the high-order corrections, and a very good 
convergence of the perturbative series at scales near = m/2. 



4. Conclusions 

Precise theoretical predictions for sparticle cross sections are essential for the in- 
terpretation of current and future searches for supersymmetry at hadron colliders. 

The inclusion of higher-order SUSY-QCD corrections reduces the scale uncertainty 
substantially. The higher-order terms also increase the size of the cross section with 
respect to the LO prediction for renormalization and factorization scales near the 
average mass of the produced SUSY particles. Thus, the SUSY-QCD corrections 
have a significant impact on the extraction of SUSY mass bounds from experimen- 
tal cross section limits, and would lead to a much more accurate determination of 
SUSY parameters like masses and couplings in the case of discovery. 

In this review we have presented the state-of-the-art SUSY-QCD predictions for 
squark and gluino hadroproduction cross sections at the Tevatron and the LHC at 
7 and 14 TeV centre-of-mass energy, including NLO corrections and NLL threshold 
resummation. We have discussed the impact of the SUSY-QCD corrections on the 
cross sections and have provided an estimate of the theoretical uncertainty due to 
scale variation, parton distribution functions, and the strong coupling as- Special 
emphasis has been placed on the predictions for the LHC at 7 TeV energy, which 
are of immediate relevance for ongoing SUSY searches. 

The effect of soft-gluon resummation is most pronounced for processes with 
initial-state gluons and final-state gluinos, which involve a large colour charge. 
Specifically, at the LHC with 7 TeV we find an increase of the cross-section pre- 
diction of up to 20% for sparticle masses around 1 TeV when going from NLO to 
NLO-I-NLL, depending in detail on the final state and the ratio of squark to gluino 
masses. Furthermore, the scale uncertainty is reduced significantly at NLO-I-NLL 
accuracy over the full range of sparticle masses relevant for hadron collider searches, 
with a remaining scale uncertainty of less than 10%. We have furthermore presented 
estimates for the overall theory uncertainty including the 68% C.L. pdf and ag error 
added in quadrature, combined linearly with the scale variation error. Even though 
the pdf and as uncertainty is significant, the inclusion of threshold resummation 
leads to a sizeable reduction of the overall theory uncertainty, which is in general 
20% or smaller for NLO-hNLL cross sections calculated with 2008 NLO MSTW 
pdfs. 

The NLO-I-NLL cross sections presented in this paper constitute the state-of-the- 
art QCD predictions for squark and gluino production in the MSSM, and provide 
therefore, we think, the optimum theoretical basis to interpret current and future 
searches for supersymmetry at the Tevatron and the LHC. 
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Fig. 1. The NLO+NLL SUSY-QCD cross section for the individual squark and gluino pair- 
production processes at the Tevatron, pp — )■ qq , qq , qg , gg + X and pp X, as a function 
of the average sparticle mass m. The error band includes the 68% C.L. pdf and ctg error, added in 
quadrature, and the error from scale variation in the range m/2 < /i < 2m added linearly to the 
combined pdf and ag uncertainty. 
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Fig. 2. The NLO+NLL SUSY-QCD cross section for inclusive squark and gluino pair-production 
at the Tevatron, pp ^ qq + qq + qg + gg + X, as a function of the average sparticle mass m. The 
error band includes the 68% C.L. pdf and Og error, added in quadrature, and the error from scale 
variation in the range m/2 < /i < 2m added linearly to the combined pdf and as uncertainty. 
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Fig. 3. The NLL /f-factor /^nll = cnlo+nll/cnlo for the individual squark and gluino pair- 
production processes at the Tevatron as a function of the average sparticle mass m. 
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Fig. 4. The scale dependence of the LO, NLO and NLO+NLL cross sections for the individual 
squark and gluino pair-production processes at the Tevatron. The squark and gluino masses have 
been set equal rriq = rrig = m in the upper four plots. 
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Fig. 5. The NLO+NLL SUSY-QCD cross section for the individual squark and jjluino pair- 
production processes at the LHC with 7 TeV, pp — > qq , qq , qg , gg + X and pp — )■ titi + X, as a 
function of the average sparticle mass m. The error band includes the 68% C.L. pdf and Os error, 
added in quadrature, and the error from scale variation in the range m/2 < fi < 2m added linearly 
to the combined pdf and as uncertainty. 
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Fig. 6. The NLO+NLL SUSY-QCD cross section for inclusive squark and gluino pair-production 
at the LHC with 7TeV, pp — > qq + qq + qg + gg + X , as a function of the average sparticlc mass m. 
The error band includes the 68% C.L. pdf and «s error, added in quadrature, and the error from 
scale variation in the range 'm/2 < < 2m added linearly to the combined pdf and as uncertainty. 
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Fig. 8. The NLL /^-factor -R'nll = ""nlo+nll/cnlo ior for inclusive squark and gluino pair- 
production at the LHC with 7TeV, pp ^ qq + qq + qg + gg + X, as a function of the average 
sparticle mass m. 
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Fig. 9. The scale dependence of the LO, NLO and NLO+NLL cross sections for the individual 
squark and gluino pair-production processes at the LHC with 7 TeV. The squark and gluino masses 
have been set equal rUq = rrig = m in the upper four plots. 
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Fig. 10. The relative NLL /f-factor i^NLL = ""nlo+nll/cnlo for the individual squark and 
gluino pair-production processes as a function of the average sparticle mass m, for various mass 
ratios r = nig/mq. 
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Fig. 11. The NLO+NLL SUSY-QCD cross section (left) and NLL ii'-factor (right) for gluino 
pair-production in the heavy squark limit, r = mg/iriq = 0, at the LHC with 7TcV, as a function 
of the gluino mass mg . 





+ qg + gg + X) 


\/5 = 7 TeV 








NLO, M=ini 
NLO + NLL, (1= li7j 
NLO + NLL, ^l = 2m 
NLO,/i = 2 m 




HQ = m 



200 400 600 800 1000 1200 
771 [GeV] 



Fig. 12. The scale dependence of the NLO and NLO+NLL SUSY-QCD cross sections for inclusive 
squark and gluino pair-production at the LHC with 7TeV, pp — > qq + qq+qg + gg + X , as a function 
of the average sparticle mass m. Shown are results for the mass ratio r = mg/mq = 1. The upper 
two curves correspond to the scale set to fi = m/2, the lower two curves to /i = 2m. 
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Fig. 13. The theoretical uncertainty for the the individual squark and gluino pair-production 
processes at the LHC with 7TeV, pp — > qq,qq,qg,gg + X, as a function of the sparticle mass 
rriq = nig = m. The error bands represent the NLO+NLL scale uncertainty in the the range 
m/2 < fJ. < 2m, and the total theory uncertainty including the 68% C.L. pdf and as error, added 
in quadrature, and the error from scale variation in the range m/2 < fJ, < 2m added linearly to 
the combined pdf and «s uncertainty. The total theory uncertainty is shown at both NLO and 
NLO+NLL. 
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Fig. 14. The NLO+NLL SUSY-QCD cross section for the individual squark and j;luino pair- 
production processes at the LHC with 14 TeV, pp qq , qq , qg , gg + X and pp — > titi + X, as a 
function of the average sparticle mass m. The error band includes the 68% C.L. pdf and Os error, 
added in quadrature, and the error from scale variation in the range m/2 < fi < 2m added linearly 
to the combined pdf and as uncertainty. 
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Fig. 15. The NLO+NLL SUSY-QCD cross section for inclusive squark and gluino pair-production 
at the LHC with 14TeV, pp — > qq + qq + qg + gg-i- X, as a function of the average sparticlc mass m. 
The error band includes the 68% C.L. pdf and «s error, added in quadrature, and the error from 
scale variation in the range m/2 < /i < 2m added linearly to the combined pdf and as uncertainty. 
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Fig. 17. The scale dependence of the LO, NLO and NLO+NLL cross sections for the individual 
squark and gluino pair-production processes at the LHC with 14TeV. The squark and gluino 
masses have been set equal rriq = mg = m in the upper four plots. 
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